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Abstract 

We consider regularized approximate cloaking for the Helmholtz equation. 
Various cloaking schemes have been recently proposed and extensively inves- 
tigated. The existing cloaking schemes in literature are (optimally) within 
I lnp|~^ in 2D and p in 3D of the perfect cloaking, where p denotes the reg- 
ularization parameter. In this work, we develop a cloaking scheme with a 
well-designed lossy layer right outside the cloaked region that can produce 
significantly enhanced near-cloaking performance. In fact, it is proved that 
the proposed cloaking scheme could (optimally) achieve p^ in M^, N > 2, 
within the perfect cloaking. It is also shown that the proposed lossy layer 
is a finite realization of a sound-hard layer. We work with general geometry 
and arbitrary cloaked contents of the proposed cloaking device. 

Resume 

Nous considerons le probleme d'invisibilite approchee pour I'equation 
d'Helmholtz. Diverses methodes ont ete recemment proposees et etudiees. 
Les techniques de quasi-invisibilite presentes dans la litterature approchent 
I'invisibilite parfaite avec une erreur proportionelle a |lnp|~^ dans et 
p dans M^, oil p designe le parametre de regularisation. Dans ce travail, 
nous developpons un systeme d'invisibilite qui utilise une couche avec perte 
a I'exterieur de la region dissimulee et ameliore considerablement la quasi- 
invisibilite. Nous prouvons que cette nouvelle technique de dissimulation 
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approche I'invisibilite parfaite avec une erreur proportionelle a p dans M , 
N > 2. Nous prouvons egalement que cette couche avec perte est un cas 
particulier d'une couche rigide. Ce travail concerne des dispositifs de dissim- 
ulation avec une geometrie generale. 

Keywords: acoustic cloaking, transformation optics, FSH lining, 
asymptotic estimates 



1. Introduction 



A region is said to be cloaked if its contents together with the cloak 
are invisible to certain measurements. From a practical viewpoint, these 
measurements are made in the exterior of the cloak. Blueprints for ma king 
objects invisible to electromagnetic waves were proposed by Pendry et al. [35 
and Leonhardt 25 in 2006. In the case of electrostatics, the same idea was 



discussed by Greenleaf et al. [16|| in 2003. The key ingredient is that optical 



parameters have transformation properties and could be pushed- forward to 
form new material parameters. The obtained materials/media are called 



transformation media. We refer to [lO|, ll3|, uA, |3J, l39|, |4J| for state-of-the-art 



surveys on the rapidly growing literature and many striking applications of 
the so-called 'transformation optics'. 

In this work, we shall be mainly concerned with the cloaking for the 
time-harmonic scalar waves governed by the Helmholtz equation. The trans- 
formation media proposed in l35l | are rather singular. This poses much 
challenge to both theoretical analysis and practical fabrication. In order 
to avoid the singular structures, several regularized approximate cloaking 



schemes are proposed in [12, [19|, |20|, |28|, |37[. The idea is either to incor- 



porate regularization into the singular transformation underlying the ideal 
cloaking, or to truncate a thin layer of the singular cloaking medium near 
the cloaking interface. Instead of the perfect invisibility, one would consider 
the 'near-invisibility' depending on a regularization parameter. Our study 



is closely related to the one introduced in j20| for approximate cloaking in 
electric impedance tomography, where the 'blow-up- a-point' transformation 
in (igI . 35 1 is regularized to be the 'blow- up-a-small- region' transformation. 
The idea was further explored in 19|, |28|, |33| for the Helmholtz equation. 
In 



28 



the author imposed a homogeneous Dirichlet boundary condition 
at the inner edge of the cloak and showed that the 'blow-up- a- small-region' 
construction gives successful near-cloak. In 19j, the authors introduced a 
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special lossy-layer between the cloaked region and the cloaking region, and 
also showed that the 'blow-up-a-small-region' construction gives successful 
near-cloak. For both cloaking constructions, it was shown that the near- 
cloaks come, respectively, within l/|lnp| in 2D and p in 3D of the perfect 
cloaking, where p is the relative size of the small region being blown-up for 
the construction and plays the role of a regularization parameter. These es- 
timates are also shown to be optimal for their constructions. More subtle 
issues of the lossy-layer cloaking construction developed in [l9| were studied 



m 
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It is worth noting that if one lets the lossy parameter in [19| go to infinity, 
this limit corresponds to the imposition of a homogeneous Dirichlet boundary 
condition at the inner edge of the cloak. On the other hand, the imposition 
of a homogeneous Dirichlet boundary condition at the inner edge of the 
cloak is equivalent to employing a sound-soft lay er right outside the cloaked 
region. In this sense, the lossy layer lining in 19| is a finite realization of the 



sound-soft lining in [28|. We would like to emphasize that employing some 
special lining is necessary for a successful near-cloaking construction, since 
otherwise it is shown in 19|] that there exists resonant inclusions which defy 
any attempt to achieve near-cloak. 

Though the existing cloaking constructions would yield successful near- 
cloaks, cloaking schemes with enhanced cloaking performances would clearly 
be of great desire and significant practical importance, especially in the 2D 
case as can be seen from our earlier discussion. A novel regularized cloaking 
scheme were developed in 2J] by making use of an FSH lining. The FSH 
lining is a special lossy layer with well-designed material parameters. The 
study in 2J] is conducted for cloaking device with spherical geometry and 
uniform cloaked contents, where the authors rely on spherical wave series 
representation of the underlying wave field to derive the estimates of the 
cloaking performance. The newly developed cloaking scheme is shown to 
produce significantly enhanced cloaking performance. In this work, we shall 
prove the general case with general geometry and arbitrary cloaked contents 
of the FSH lining construction. For the construction, it is shown that one 
could achieve, respectively, p^ in 2D and in 3D within the perfect cloaking. 
Apparently, our novel cloaking proposal with such significantly improved 
cloaking performances would be a very promising scheme for constructing 
practical cloaking device. From our arguments in deriving these estimates, 
one can see that the FSH layer is a finite realization of a sound-hard layer. 
Hence, the FSH layer is of completely different physical nature from the one 



3 



m 



19| which is a finite reahzation of a sound-soft layer. In fact, the one in 



19[ makes essential use of a large lossy parameter, whereas for our FSH layer 



we only require a finite lossy parameter but a large density parameter of the 
layer medium. 

Th e analysis of cloaking must specify the type of exterior measurements. 

In 



12, 19, 20 



the near-cloaks are assessed in terms of boundary measure- 
ments encoded into the boundary Neumann-to-Dirichlet map or Dirichlet-to- 
Neumann (DtN) map. The scattering measurement encoded into the scat- 



tering amplitude is considered for the near-cloaks in [24, |28| . In the current 



article, we shall assess our near-cloak construction with respect to the bound- 



ary measurements. Nonetheless, by |40l . |41| . it is known that knowing the 



boundary DtN/NtD map amounts to knowing the scattering amplitude. 

In this paper, we focus entirely on transformation-optics-approach in con- 
structing cloaking devices. But we would like to mention in passing the other 
promising cloaking schemes including the one based on anomalous localized 
resonance [3l|, and another one based on special (object-dependent) coatings 
[l|. It is also interesting to note a recent work in where the authors im- 
plement mult i- coatings to enhance the near-cloak in EIT. The same idea has 
also been extended to acoustic cloaking for achieving enhancement in 0, [sj. 

The rest of the paper is organized as follows. In Section 2, we develop the 
cloaking scheme by employing the FSH lining and present the main theorems. 
Section 3 is devoted to the proofs of the main results. In Section 4, we derive 
some crucial estimates on small inclusions that were needed in Section 3. In 
Section 5, we consider our cloaking construction within spherical geometry 
and uniform cloaked contents, which illustrates the sharpness of our estimates 
in Section 3. Section 6 is devoted to discussion. 



2. Near-cloak with FSH lining 

Let q G L°°(M^) be a real scalar function and a = {(7^^)fj^^ G Sym(A^) 
be a symmetric-matrix- valued function on M^, which is bounded in the sense 
that, for some constants < cq < Co < cxd, 

coei<e^{x)i<c,ei (2.1) 

for all X G and ^ G M^. In acoustics, and g, respectively, represent 
the mass density tensor and the bulk modulus of a regular acoustic medium. 
We shall denote {M^; cr, q} an acoustic medium as described above. It is 
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assumed that the inhomogeneity of the medium is compactly supported, 
namely, a = I and g = 1 in ]R^\Z) with D a bounded Lipschitz domain in 
M^. In M^, the scalar wave propagation is govern by 



qix)Uu - V • (or(x)Vf/) = in 



The time-harmonic solutions U{x,t) = u{x)e is described by the hetero- 
geneous Helmholtz equation 

i,j=l * 

Let C be a bounded Lipschitz domain such that D C Q. An important 
problem arising from practical applications is described as following. Let 
^ = {^i)iLi be the exterior unit normal vector to dfl. Impose the following 
boundary condition on dQ for (12. 2p . 

Y u,a'^ ^ = ^ e H-^/\dQ) on dQ, (2.3) 

and define the Neumann-to-Dirichlet (NtD) map by 

K^,{i;)=u\9n^H^'\d9), (2.4) 

where u G H^{VL) solves (E2])-(E3]). It is known that A„,g : H-^/^{dQ) 
H^/'^{dVL) is well-defined and invertible provided avoids a discrete set of 
eigenvalues. The practical problem is to recover {D; a, g} by knowledge of 
Ao-^g which encodes the exterior boundary measurements. 

In this paper, we shall be concerned with the construction of a layer 
of cloaking medium which makes the inside medium invisible to exterior 
measurements. To that end, we present a quick discussion on transforma- 
tion acoustics. Let x = F{x) : Q ^ Q he a. bi-Lipschitz and orientation- 
preserving mapping. For an acoustic medium {Q;a,q}, we let the push- 
forwarded medium be defined by 

{fi; a,q} = F,{n; a, q} := {Q; F,a, F,q}, (2.5) 

where 

a{x) = F,cr{x) := jMcr{x)M^\,=F-^i) 
q{x) = F^q{x) := q{x) / J\^=f-\x) 
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and M = {dxi/dxj)f^j^^, J = det(M). Then u E solves the Helmholtz 

equation 

V ■ (o"(x) Vu) + u'^q{x)u = on fi, 
if and only if the pull-back field u = {F^^)*u := u o F^^ G H^{Q) solves 

V ■ {a{x)Vu) + uj'^q{x)u = 0. 

We have made use of V and V to distinguish the differentiations respectively 
in x- and x-coordinates. We refer to |19|, |28| for a proof of this invariance. 

We are in a position to construct the cloaking device. In the sequel, we 
let Qhe a. connected smooth domain and D be a convex smooth domain, and 
suppose that D <^ Q and Q\D is connected. W.L.O.G., we assume that D 
contains the origin. Let p > be sufficiently small and Dp := {px;x G D}. 
Suppose 

Fp : n\Dp n\D, (2.7) 

which is a bi-Lipschitz and orientation-preserving mapping, and Fp\gQ = 
Identity. A celebrated example of such blow-up mapping is given by 



y = Fp{x) := ( -^^R2 + -i, -\x\ ) — , p < Ri < R2 (2.8) 



' Ri- p R2- Ri. \ X 

-B R2 + -B fI y- 

R2 - p R2- P / F 

which blows-up the central ball Bp to Br^ within Bji^- Now, we set 



, , ■ Fp(x) for X G Vl\Dp, , ^ 

F[x) = { ^ , ^ " (2.9) 



for X G Dp. 



Clearly, F : i7 — > f2 is bi-Lipschitz and orientation- preserving and F\qq^ = 
Identity. Next, let 

{Dp\Dp/2;cri,qi}, = jp^+'l, qi = a + i(5, (2.10) 

where a, /3, 7, 5 are fixed positive constants, and 

{D\D,/2-, o[, i} = F,{Dp\Dp/2; ai, qi}. (2.11) 

We further let 

{n\D; <, q^J = iFpUn\Dp; /, 1}. (2.12) 
Let -Di/2 represent the region which we intend to cloak and 

{Dy2;K,q'a} (2.13) 
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be the target medium which is arbitrary but regular. We claim the following 
construction yields a near-cloaking device occupying Q, 



aP,qP mn\D, 
{n;a,q}={ a'i,q[ m D\Dy2, (2-14) 
Qa in ^1/2- 

In order to present the main theorem justifying the near-cloaking con- 
struction f l2.14p . we let uq be a solution to the following PDE system 

Auq + u^Uq = in 

duo , (2-15) 
= w on o\l. 

ov 

That is, Uq is the wave field in the "free space". We suppose that —oj"^ is 
not an eigenvalue of the Neumann Laplacian. Hence, one has a well-defined 
"free" NtD map 

Ao(V') = Molan, 
where Uq solves fl2.15p . We have 

Theorem 2.1. Suppose —bJ^ is not an eigenvalue of the Laplacian on Q with 
Neumann boundary condition. Let Ao-,^ be the NtD map corresponding to the 
construction ^2.14\ ), and Aq be the "free" NtD map. Then there exists a 
constant po such that for any p < po, 

||A^_g - Ao||^(^-i/2(gQ)^^i/2(gQ)) < Cp^, (2.16) 
where C is a positive constant dependent only on po, u, a, (3, 7 and D, VL, 



but completely independent of p. That is, the construction ( 2.14 ) produce a 



near-cloaking scheme which is within p^ of the perfect cloaking in M^. 

3. Proof of the main result 

This section is devoted to the proof of Theorem 12.11 First, for {Q]a,q} 
given in (12.141) . we let 

in Q,\Dp, 

{Q;(rp,qp} = {F'^)^{n](r,q} = \ ai,qi mDp\Dp/2, (3.1) 

in Dp/2, 
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where 

{Dp,^-<Ja,q,} = {F-%{D„^-a'^,q'^}. 



We consider the solution Up of 



V ■ {(JpVup) + oj'^qpUp = in fi, 

, (3-2) 
— %p oil. 



dv 

Noting F\q£i = Identity, by the transformation acoustics, it is straightforward 
to show that 

A.,,(^)=A.^,,^(^), y^p e H-'/\dn). (3.3) 
Hence, in order to prove Theorem 12.11 we only need to show 

Theorem 3.1. Suppose —u"^ is not an eigenvalue of the Laplacian on Q with 
Neumann boundary condition. Let uq and Up he the solutions of Ii2.15\) and 
lis. 2^) respectively. Then there exists a constant pq > such that for any 
P < Po, 

hp - Mo 1 1 J/1/2 (af7) < Cp^\\ilj\\H-i/2^an), (3.4) 

where C is a constant dependent only on po? ^> oi, 7 ^^^'^ D, Q, but 
independent of p and tp. 

Our proof of Theorem 13.11 would follow the spirit of the one for proving 



the main theorem in [19[. However, the main strategy in [19| is to control 
the Dirichlet value of Up on the exterior of the lossy layer, namely dDp, and 
then derive some estimates of exterior boundary effects due to small sound- 
soft like inclusions; whereas in our case, we would control the value of the 
conormal derivative of Up on the exterior of the lossy layer dD'^, and then 
derive some estimates of exterior boundary effects due to small sound-hard 
like inclusions. It is also emphasized that by making use of layer potential 
techniques, we work with general geometry of the cloaking device. 
We first derive the following lemma. 

Lemma 3.2. The solutions of i\2.15\) and liS. 2\) satisfy 



I3uj _ \up\ dx < C\\^Ij\\h-i/2(^q^)\\up- uoWH^^f^an)^ (3-5) 
where C is a positive constant (depending only on Q). 
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Proof. Multiplying fl3.2p by Up and integrating by parts, we have 

— / ap\Vup\^ dx + uj^ I qp\up\^ dx = — (cTpVup) ■ uUp da^, 
Jn Jn Jon 



(3.6) 



which in turn yields 



f3uj / \up\ dx 

Jd2p\Dp 

= ( [ ^ ■ Up da^ ] = [ I il){up - Mo) da. 



(3.7) 



By (13. 7p . we immediately have (13. 5p . 
In the following, we let 

m-{x) = u-Vu-p{x) on dDp, (3.8) 

namely, the normal derivative of Up on dDp when one approaches dDp from 
the interior of Dp. Here and throughout the rest of this paper, u denotes the 
exterior unit normal of the domain under discussion. Similarly, we let 

^+(x) = 1^ ■ Vu^{x) on dDp (3.9) 

denote the normal derivative of Up on dDp when one approaches dDp from 
the exterior of Dp. We shall show 

Lemma 3.3. The solutions to Ii2.15\) and liS. 2\) verify 



3/2(9D) 



^ Al±V^±W^^ -N-2n , II II II ^^-^^^ 

< ^^2^2 P lFllH-i/2(on)lFp - wo||Hi/2(an), 



and 

||^+(P ■)||i-3/2(a^) 

^ — P lFilH-i/2(9n)||Mp - Mo||//i/2(an), 

where C is positive constant dependent only on D and Q, but independent of 
ip and p. 
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Proof. We shall make use of the following fact 



W^^ip ■)\\H-3/HdD)= sup / ^{px)<j){x) da, . (3.12) 



f3/2{SI5)-^ 

For any G H^^'^{dD), there exists w E H'^(D) such that 
(i) w = (j) on dD and ^ = on dD, 



(ii) \\w\\h^D) < C||0||^3/2(g^), 



[iii) w = in D 



'l/2- 



Then we have 



len we nave 

f f du- f du 

/ ^ (px)(/)(x) da^ = / ——{px)(t){x) d(T,= — 

JdD JdD JdD 



X := - e D. 

P 



v{x) := Up{px) = Up{y), x E D. 



For y e Dp, let 

Set 
Since 

'fVy ip'^^^VyUp)+u\a + il3)up = in Dp\Dp/2, 
it is directly verified that 

-fV,- (p^V^v) +uj^{a + i(3)v = in D\Di/2. 

By Green's formula and f l3.13p . we have 



{px)w{x) da,. 

dD 

(3.13) 



(3.14) 
(3.15) 



\1' {px)(f){x) da. 



dD 



— ^(px)0(x) da,= p ^ 
dD ov Jqd dv 

1 



9f " 



(X)</)(X) rfdj; 



■{x)w{x) da, 



/ ^(a;)^(a;) rfcxa; 
dD du 



(3.16) 



,-1 



Af (x)w(x) dx — v{x)Aw{x) dx 

D Jd 
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Then by fl3.15p and fl3.16p . we further have 



^ {px)(f){x) da^ 



dD 



< 



Av(x)w(x) dx — v(x)Aw(x) dx 

D J D 

V^^TW -S-l 2[ f I , ,|2 , V^'ll II 

P w / \v{x)\ dx \\w\\L2m) 

+ p-M / _ \v{x)\^dx] \\Aw\\L2^Dy 

Using the relation 

\MlHd\d,/,) = \Mp ■)\\lHd\d,^,) = P'^'^'l' 
we have from (13.171) that 



\1/ (px)0(x) da^ 

'dD 



(3.17) 



u 



■p\\L\Dp\D^^2) 



(3.18) 



'<5, .2 



which imphes 

P'ip ■)L-3/2(9D) 
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< Cp- 



N/2-1 ( 1 ^ y^^+Ep-sJ 



\U 
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pllL2(Dp\D^/2)- 

(3.19) 

By (I3.19P and Lemma 13. 2[ one immediately has (I3.10p . Finally, by (I3.14p 
and the transmission condition on dD, we see 



du~ 



du 



IP 



2+5 



du 



dDp 



dp 



dDp 



and hence 

vl/+(px) = 7p^+'^^"(pa;) for x & dD 

which together with (I3.10p implies (13. lip . 
The proof is completed. 
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The next lemma is of crucial importance in proving Theorem 13.11 

Lemma 3.4. Suppose —bj^ is not an eigenvalue of the Laplacian on Q with 
Neumann boundary condition. Let uq € H^{Q) be the solution of Ii2.15\) . Let 
E H^^/'^[dDr) and consider the Helmholtz system 



Am^ + uj'^Ur = in ^l\Dj 

An 

(f on oDt, 



V 



(3.20) 



dp 



Let 



dup 
du 



on dQ. 



[x) for X G dDj 



Then there exist a constant tq > such that for any r < tq, 



H-i/2(aQ) 



+ r 



N-l 



-nH-y^aD)) , (3.21) 



where C is a positive constant dependent only on Tq, u and Q, D, but inde- 
pendent of T and if, ip. 



Proof. Let 

V = Ur — uq on VL\Dr. 
By fl27[5D and ([3201), one sees that V G H\Q\Dr) satisfies 



(3.22) 



AV + u^V = in n\D, 

dV duo ^„ 

— = Lf- — — on dD^, 

du 

— = on dn. 
du 



Let V = Vi-V2 with Vi G H\n\Dr) satisfying 



(3.23) 



AVi + cu^Vi = in n\D^, 



dVi 
du 







on dDr, 
on dil. 



(3.24) 
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and V2 E H^{Q,\Dr) satisfying 



AV2 + = in n\Dr, 
dV2 



du 







on dD^-, 
on dil. 



(3.25) 



By Lemma 14.11 in Section HJ we know 



l^2||/fi/2(an) < Cr^ 



/i'-i/2(an)- 
1 mAfv 



(3.26) 



In order to estimate ||K||hi/2(oq), we let W G Hl^^{M. \Dt-) be the unique 
solution to the following scattering problem 



iioW )■ = 0. 



AW + u^W = 
dW 

hm 



Let To be sufficiently small such that 



in R^\D 
on 



(3.27) 



(3.28) 



for some finite tq > 0, where Bj. denotes a central ball of radius r. Let 
ro < ri < r2 < +00 be such that 



Br, m n and n\D d Br^\Br,. 
By Lemma 14.21 in Section HI we have 

Since (A + u'^)W = 0, by the interior regularity estimates, we see 



(3.29) 



(3.30) 



dW 




dv 





<Cr^^-iv^(r-) 11^-3/2(5,,), 



(3.31) 
(3.32) 
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and 



l^llHi/2(af7) 



Next, by the Green's representation, we know 



(3.33) 



W{x) 
where 



^-^^^W{y)-G{x-yf-^day, x G (3.34) 



G{x) 



, \ (^-2)/2 



4 V27r x 



H. 



(7V-2)/2 



'uj\x\) 



(3.35) 



_ _dW 

du 



dv 



By Lemma 14.31 in the following, we have 



\P\\m/-^{dn) - ^ 



dW 




du 


an 



which together with f l3.36p implies 

Since Vi = P + W, fl333|) and f l339|) immediately yields that 



which together with fl3.26p implies fl3.2ip . 
The proof is completed. 



(3.36) 



is the outgoing Green's function. By (I3.3ip . fl3.32p and fl3.34p . it is readily 
seen that 

^'^ an c{an) 

Let 

p = Vi-W. 

By (Km and (Km . one sees that P G H\Q\Dr) satisfies 

AP + u'^P = in n\Dr, 
on (9-Dt-, 

on d^l. 



(3.37) 



(3.38) 



(3.39) 



14 



We are in a position to present the proof of Theorem I3.1[ 



dut I 



Proof (Proof of Theorem 13. II ). By taking r = p and ip = -g^laop in 
Lemma [3.41 we have 



(3.40) 



Next, by (13. lip in Lemma [3.31 we have for e > 



du 



(P-; 



/f-3/2(eD) 



1 1/2 



<C2P 
+ 



I J2-N)/2 . N-1 
^(2-N)/2{ P P 

4e 



(3.41) 



P' 



{2-N)/2 . pN-1 II^P ~ ^0||//i/2(9n) 



From f l3.40p and (13.41 p . we further have 

hp - %||//i/2(0n) <C'lP^II^II//-i/2(9Q) 



+ ^C'lCap^ ^^ll^ll//-i/2(0Q) + CiCaellwp - Mo||iTi/2(ao). 

(3.42) 

By choosing e such that CiC2e < 1, we see immediately from (I3.42p that 

hp - %||Hi/2(9n) < C'P^ll^llH-i/2(an)> 
which completes the proof. 



4. Some estimates on small inclusions 

In this section, we shall derive those lemmas that were needed in the proof 
of Lemma 13.41 on the wave estimates due to small inclusions. We would like 
to mention that there are a lot of results on this subject in different settings 
in literature, see e.g., 0,0,0 , [oj . We shall derive some new estimates in the 
specific setting of our current study. We would make essential use of the layer 
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potential techniques to derive the desired estimates in this section. To that 
end, we let G{x) be the outgoing Green's function in (13.351) . It is well-known 
that when N = 2, 

G(x) = -— In Ixl + - - — In - - — + In Ixl) (4.1) 

^ ' 27r ' ' 4 27r 2 27r VI I I i; v ; 

for |x| — 0, where E is the Euler's constant; and when = 3 

G{x) = ^ (4.2) 

For surface densities iplx) with x G dQ, and (f{x) with x G dD-r, we introduce 
the single- and double-layer potential operators as follows 



and 



{SL[9n]^){x) = [ G{x- y)ij{y) day, x G M^\5fi 
Jan 

{SL[9D^]^){x) = [ G{x- vMy) day, x G 

JdDr 

{DL^en]i^){x)= [ ^^j'^'/V lz/) day, x G M^\an 
Jan (J^KU) 

{DL^9DAV^){x) = j^^ ^^^^^{y)day, xGR^\9D.. 



(4.3) 



(4.4) 



We also let 



and 



{S[an]'4^){x) = / G{x - y)ip{y) day, x e dQ 
Jan 

{S[OD^]^){x) = / G{x - y)ip{y) day, x G dD^ 
Jan^ 

r dGix - y) 
{K^aa]^){x) = ^^^^^ — ip{y) day, x G dVL 

r dGix - y) 
{K^QD,]'^){x) = / — ip{y) day, x G dD^. 



(4.5) 



(4.6) 



I an. My) 

As specified earlier in Section 3, we would like to emphasize again here 
that in the above integral operators, u denotes the exterior unit normal vector 
of the underlying domain. 
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Lemma 4.1. Suppose —iJ^ is not an eigenvalue of the Laplacian on Q with 
Neumann boundary condition. Let uq G H^{n) be the solution of ^2.15\) and 
(po{x) = ^^{x) for X G dDr- Consider the Helmholtz system 

Aur + uj'^Ur = in 0,\Dr, 
— = ipo on dDr, (4.7) 

— — = on oil. 

Then there exists a constant tq > such that for any t < tq, ( [^. ?] ) has a 
unique solution Ur G H^{Q\Dr) and moreover 

\Mm/2(^9n) < C II II 1/2 (an)- (4-8) 

where C is a positive constant dependent only on tq, u and Q, D, but inde- 
pendent of T and ip. 

Proof. Since (^n G C{dDr), we know G C^(fi\D^) nC(fi\D^) is a strong 
solution (cf. By Green's representation formula, we know 

Ur{x) = [ \g{x -y) ^!!''!^^} - ^^i^. J^\ r{y)\ day, Xen\Dr. 

Jd{n\Dr) I duiv) duiy) J 

(4.9) 



Let 



and 



h{x) := - / G{x - y)(po{y) dcTy 

'dDr 



01 = Ur\dD^ and 02 = Ur\dn- 

From (14. 9 p we have 

u,{x) = {DL[9D^]<f)i){x) - {DL[9n]<j)2){x) + h{x), x G Q\Dr. (4.10) 

Letting x go to dQ and OD^-, respectively, by the mapping properties of 
double-layer potential operator (cf. jlH and jsoj), we have from (I4.10p the 
following system of integral equations for 0i G C{dDr) and 02 G C{dQ), 

\(pi{x) = (fCrai) i0i)(x) - {DLi9ai(f)2){x) + h{x), x G dDr 

(4.11) 

-02(x) = {DLi9Dr](pi)ix) - {K[9Q](j)2)ix) + h{x), X G dn. 
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Next, we claim 



(4.12) 

We first estimate h{x) for x E dfl. It is noted tliat \G{x — y)\, \VyG{x — y)\ < 
C for any x G dQ and y G dDr, where C depends only on tq and dQ. 
Integrating by parts, we have for any x G dQ 



\h{x)\ 



^/ .duoiy) 



/ AyUo{y)G{x - y) dy + VyG{x - y) ■ VyUo{y) dy 

J D-r J D-r 

/ S/yG{x - y) ■ VyUoiy) dy 

J 



G{x - y)uo{y) dy 

<Cr^ (w^||^io|U°°(DO + ||Vmo||l°=(do) 
<Cr^||^||^-v.(an)- 



(4.13) 



We proceed to estimate h{x) for x G 00^. By taking tq sufficiently small, we 
assume D2ro fi. We first let x G SD^^ with 1 < 5 < 2. Similar to fl4.13p . 
by integration by parts, we have 



\h{x)\ < 



G{x - y)uo{y) dy 



Dr 



/ VyG{x -y) -VyUQiy) dy 

J Dr 



For the first term in (14.141) . we have for x' G dD^ 

\h,{Tx')\<uj^ I \G{t{x' - y'))u,{Ty')\T'^ dy' 



(4.14) 



(4.15) 



D 



Using (14.11) and (14.21) . and the mapping property of volume potential opera- 
tor, it is straightforward to verify that 



\hi{T ■)\\c{dDs) < Ct\\uo{t OIIl-cd) < C^||^||H-i/2(an), 



(4.16) 



where G is independent of 6 and r. In like manner, it can be shown that for 
the second term in (I4.14p 



h-2(x) 



VyG{x - y) ■ VyUo{y) dy, 



Dr 
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we have 

\\h2{r ■)\\cidD,) < CrllV'llH-i/^can)- (4-17) 

By the mapping properties of single- layer potential operator (cf. [11]), we 
know 

h{x)\aD, = lim {h{x)\aDs^) , 

5— s>l+ 

which together with fl4.16p and fl4.17p implies 

\\h{T-)\\ci9D)<CT\mH-i/2^9^y (4.18) 

We now consider the system of integral equations fl4.1ip . It is first noted 
that 

II {DL^anM (■)\\cion) < Cr^'^-'^/'UM-nLHaD.)- (4-19) 

It is further noted that since —co^ is not an eigenvalue of the Laplacian on Q 
with Neumann boundary condition, 

h + K[9n] : c{dn) ^ cidn) 

is invertible (cf. [ll|). Hence, by the second equation in (14. lip , and (I4.12p . 
(I4.19p . we have 

||02||c(an)<C'(r(^-i)/2||0i(-)IU^(aD.)+r^||^||H-v.(ac)). (4.20) 

We proceed to treat the first equation in ( 14. lip . First, by change of variables 
in integrals, it is straightforward to show that 

{K[9D^]<Pi) {rx') = {KoiaoMr •)) (™') + Wi{r ■)) (™'), e dD, 

(4.21) 

where -R'opD] is an integral operator with the kernel given by dGo{x' — y') / dv{y') 
with 

ln|a;'-y'| = 2, 

1 (4.22) 

r A^ = 3, 

47r \x' — y'\ 

for x' 7^ y' and x' ,y' G dD, and satisfies 

ciL^idD),L^iaD)) < C'e(r), (4.23) 
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where 

J In r when N = 2, 
~ [t"^ wheniV = 3. 

Using fl4.2ip . the first equation in (14.1 ip can be reformulated into 
1 



irx') + {DL[9n]M-)) (™') = K^x'), x' e dD. 

(4.24) 

Then using (I4.20p . it is directly verified that 

Since /- iKo[aD] is invertible from L'^{dD) to L'^{dD) (see |42j, [381, §7.11]), 
by (OS]) . dOl]), (I4:25|) and (HA2|) . one can show that 

\\Mr-)\\LH0D)<CT\mH..^.^9^y (4.26) 

Using the relation ||</)i(-)IU2(aD,) = r(^"^)/2||0(r ^[[^^(az)), one further has 
from fICTll that 

||</'i||L^(aD.)<C^^^'''^/'ll^ll/z-/.(an). (4.27) 
Then, by Km and (S2ID, we see that 

||02||c(an) < C^r^||^||H-V2(en)- (4-28) 
Finally, by the second equation in (14. lip , we have 

uix) = 2 {{DL[^D.](|)^){x) - iK[9n](l>2)ix) + h{x)} , x E dn. (4.29) 

By a similar argument to (I4.13p . one can show that H/iUcH^f^) — C'^^\\'^\\H-'^/^{dQ)j 
which implies 

\\h\\HyHon)<Cr''\mH-yHdny (4-30) 

Since for x G dfl, DLi9d-^] has a smooth kernel, by (14.27^ and Schwartz 
inequality, it is straightforward to show that 



\DLi9d^](I)i\\hi/2(^9q) < C[mecisuYe{dDr)Y^^\\(t)i\\L2(9D^) < Cr^\i(fj\\H~i/2(^9n)- 

(4.31) 
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By noting that K]^qq^ is bounded from to H^{dQ) (cf. [ll|), we see 

from IK^ that 

\\Kian](l)2\\m/Hdn) < CT^Mn-^/Han)- (4-32) 
Combining f lOg]) -f O^ . we have 

ll'"ll//i/2(aQ) < Ct^\\iIj\\h-i/2^9q-), 
which completes the proof. 

Lemma 4.2. Let B^J = 0, 1, 2, D, be the ones in i\3.28\) and \3.29\) and 

let W G HI^^(R^\Dt-) be the unique solution to 

AW + u^W = m R^\Dr, 

-^ = ^eH-'/'idD^) ondDr, (4.33) 

|a;|-5-+oo lo\x\ J 

There exists a constant tq > such that for any t < tq, 

\\WhHB.^\B.,) < CT^'-'UiT OIIh^/.^d), (4.34) 

where C is a positive constant dependent only on D and tq, uj, ro,r2, but 
independent of (p and r. 

In order to gain more insights, we first present a proof of Lemma 14.21 
within spherical geometry. That is, we shall first assume that D^^ = B^, the 
central ball of radius r in M^, and D = Bi. We shall make essential use of 
series representation of the wave field W. 

Proof (Proof of Lemma [Ol) . Clearly, in order to show (14. 34 p . it suf- 
fices to prove that there exist two constants Ci, C2, such that for all 0(r-) G 

H-^'^{dBi) and r < Tq, 

W\\m/\dBr^) < C'2r^"l0(r-)||//-3/2(aBi). 

We shall make use of eigenvalues and eigenfunctions of the Laplace- 
Beltrami operator on a sphere dBj. to define the Sobolev space H^{dBj.) 
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for r > and s G M, which we briefly review in the following and we refer 



32I . §5.4] and [26|, §1.7] for general discussions. The Laplace-Beltrami 



to 

operator on a circle dB^ in two dimensions is 

1 d^u 



r2 



where (r, 9) is the polar coordinate in (^q{_ page 234]); and the Laplace- 
Beltrami operator on a sphere dB,. in three dimensions is 



^dBrU 



sin 6 



d , . ^du. 1 d'^u 
oO o9 sm 9 0(f'' 



where (r, 6', (y?) is the spherical coordinates in (see |36|, Appendix] and 27|). 
By direct calculations, one has 



An9 



-A, 



^2 gine 



dBr 



V27fr r2 727rf 



2 



/2^r}n^-oo is an orthonormalized basis in L {dBr) and A„ := /r 



is the eigenvalue corresponding to the eigenfunction w„. Suppose u{x) 
Er=-ooCn(r)e*"^ e H'{dBr), s G M, then by we know 



H^idBr) 



E (1 + 



)y27rr| 



(4.36) 



It is noted that when r = 1, this is consistent with the if'^[0,27r] presented 
in 21I, §8.1]. In three dimensions, by straightforward calculations, we have 



-A 



K™(x) 



dBr 



-^n(n + 1 



where Y^{x) for x G n G NU{0} and m = —n, — (n — 1), . . . , {n — l),n, are 
the spherical harmonics. Hence in 3D, the eigenvalues and eigenf unctions are. 



respectively, A„ = n{n + l)/r and 



m 



-n, -(n - 1), (n 



l),n. Suppose u{x) = En=oEl=-n<ir)Y;riS:) G i/^(aS,), s G E, then 
we have 

nin + 1 



n=0 m,=— n 



(4.37) 
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We first consider the two-dimensional case of tlie lemma. Suppose r < 
To < ''^0/4 and (j){x) = (pijx') = X]^-oo '^"(''")^™^' ^ ~ tx',x' G dBi. Since 
0(rx') G H-^/^{dBi), one first sees from (ICTj) 



Suppose that the solution is of the form W{x) = J2'^=-oo^'n-H^\^\^\)^ 
X G M.'^\Br. According to the PDE system (14.33^ . we have a„ 
and 



And 



n=—o 

Hence, by direct calculations 



n=— 00 



/2 



v27rro 



< E (1 + ^')"'^'27r|c„p (1 + 5^ (1 + 



2\2 



(4.38) 



where we made use of the fact that (1 + /rl) < (1 + l/rQ)(l + n^). By 
f l4.37p and fl4.38p . it is sufficient for us to study the asymptotic development 

of 



T{r) ■.= {l + l/riy/' (1 



n 



2\2 



Since H^l{ujr) = {—l)^Hn\oor), we only need consider n > in the above 
series. By the asymptotic behaviors of Hankel functions as r — )■ +0 or 
n — )• +00 (see ^ and [i^), it can be shown that there exists a positive 
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integer Nq such that 



r, n = 0, 



(4.39) 



T T 
-(- 

n ro 



n> N< 



0- 



Furthermore, by our assumption on r, we have 

-,2\2 / _ \ 2n 



n>No 



n>No 



42r, 



(4.40) 



Using (14.391) and (14.401) . one can show by straightforward calculations that 
T{t) < Cfr^, where Ci is independent of r for tq sufficiently small. There- 
fore, we have 

< C^r\\<P{r■)\\ 

H-3/2{dBi)- (4-41) 

In like manner, one can show that there exists C2 such that 



(4.42) 



It is interesting to remark that by similar arguments, one can actually show 
that for any positive integer s, we have 



and 



m/^dBr^) < C2T\\(f){T-)\\H-s/2^QB,y 



(4.43) 



(4.44) 



where Ci and C2 are positive constants dependent only on s, tq, r2, tq and 
u, but independent of r and (f). 

The three dimensional case can be proved similarly. Suppose 0(x) e 
H-^/\dBr), and 

00 n 
n=0 m=— n 
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Then, 



oo n 



n=0 m=—n 

Suppose the solution is of the form 



oo n 



n=0 m=—n 



-. Let 



By using the boundary condition on dB^-, we have =~ 

dn = {1 + n{n + 1)) in the following calculations. It is noted that (1 + n{n + 
l)/^o) < (1 + ^/rl)dn. We have 



oo n 
n=0 m=—n 



< 



. n=0 m=~n 



3/2, rn. M2 




n=0 



(1)' 

/inHwro)ro 



(4.45) 



By the asymptotic properties of the spherical Hankel functions hn^ (lot) and 
their derivatives as r — )■ +0 or n — )• +oo (see ^ and j^sj), similar to the two 
dimensional case, one can show that 



n=0 



hn\u}ro)ro 



(1)' 



< CtT\ 



Therefore, we have 

\\W\\m/2(^9Br,) < CiT^UiT-)\\H-3/^dB,)- 

In like manner, one can show 

1 1 11^1/2(05^^) < C2T^\\(l){r-)\\H-3/2^QBi)- 

Similar to the estimates in fl4.43p and ( I4.44p for the two dimensional case, 
one can derive more general estimates for the three dimensional well. 

Next, we shall present the proof of Lemma [4.21 within general geometry, 
which is based on layer-potential techniques. 
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Proof (Proof of Lemma [¥7^ . Let 

= W\aD. G H^'\dD,). 
Similar to firTOj) . we have (cf. ^]) 

= (DL[azp.]<^)(x) - (5L[aB^]</>)(x), a; G R^\D.. (4.46) 



By the jump properties of layer potential operators (cf. [30|), we have from 
dOH]) that 

\p{x) = {K[QD^^^){x) - (S[aD.]0)(a;), x G dD^. (4.47) 

Next, we only consider the 3D case, and the 2D case could be shown in a 
similar manner. 

We shall first show that 

Mr -nH-^moD) < CrUir OIIh-^/^c^d)- (4.48) 

To that end, we first note that by straightforward calculations 

iS[9D.^(P){Tx') = T{So[aD]^{T ■)){Tx') + {y^{T ■)){tx'), x'edD, (4.49) 

where 5*0 [^d] is an integral operator with the kernel given by Go{x' — y') as 
the one in (14.221) . and ^ is an integral operator with the kernel given by 



iu 2 I T{iujTY \x' — y'\ ^ T{iujT)^ \x' — y''^ 



^(x'-y') = — r , , , , 

^ ^ ^ Air 2! 47r 3! An 



-A{\x'~y'\), (4.50) 



where A{t) is an (real) analytic function in t G M. Hence, by the mapping 
properties presented in [ii', §4.3], one has 

||(^0(r-))(ra;')|| 
Moreover, we know (cf. (ssl, §7.11], jsil, §4.4]) 

<C||0(r Oil H-3/2{dD)- (4-52) 

Next, by using the decomposition (14.211) of K^^qd^^, (14.471) can be reformu- 
lated as 



[TX 



(4.53) 



--r{So[9D]^{T ■)){tx') + {y<P{T ■)){tx'), x' G dD. 

26 



By straightforward asymptotic expansions and also using the mapping prop- 
erties presented in [32|, §4.3], one can readily show that 



£(//-i/2(eD),H-i/2(aD)) 



<Ct\ 



(4.54) 



We shall also make use of the fact that (cf. jssi, §7.11]) 



1/ 

2 



is an isomorphism from H ^^"^(dD) to H ^^^(dD). (4.55) 



Hence, by combining fl4.49p - fl4.55p . one readily has (14. 48 p . 

Finally, by taking x e -BraV-Bro in f!4.46p and using fl4.48p . we have fl4.34p 
by straightforward verification. 



Lemma 4.3. Let D,Dj- and n be the ones in Lemma 4-1 ond let P G 

H\n\Dr) satisfy 



AP + uj^P = 
dP 







in fi\-Dr5 
on dDr, 



P 



dv 



(4.56) 



^ e C{dn) on dn. 



Suppose —iJ^ is not an eigenvalue of the Laplacian on Q with Neumann 
boundary condition. Then there exists a constant tq > such that for t < tq, 
( 1^.5^ is uniquely solvable and satisfies 



l-P||//i/2(af7) < C||v2||c{9n), 



(4.57) 



where C is a positive constant dependent only on Tq, u and Q, D, but inde- 
pendent of T and (p. 

Proof. We shall make use of layer potential techniques again to show that 
lemma. The argument would follow a similar spirit to that for proving 
Lemma 14.11 but would be comparatively simpler and we shall only sketch 
in the following. 

Clearly, P e C^{^l\Dr) n C{Il\D^) is a strong solution. By letting 



we have 



Pi = P\dD^ and p2 = P\dn, 
P{x) = {DL[9Dr]Pi)ix) - iDL[gn]P2){x) + g{x) 



(4.58) 
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where 

g(x) = / G{x - y)ip{y) day 
Jon 

satisfying 

\\g{r ■)\\c(dD) <C\\(f\\c{dn) and \\g\\c{dn) < C\\v\\cidn)- (4.59) 

By the jump properties of double-layer potential operator, we have from 
fl4.58p the following system of integral equations for pi G C^dD^-) and p2 G 
C{dn), 

lrpi{x) = {KiQD,]pi){x) - {DLion]P2)ix) + g{x), x G dDr 
I (4.60) 
-P2ix) = {DL[9D^]pi){x) - {K[9n]P2)ix) +gix), xedVt. 

By a similar scaling and asymptotic argument to that in the proof of Lemma I^T 
one can show that 

\\Pi\\c{dDr) <C\\Lp\\c{dQ.) and \\p2\\c{aQ.) < C\\'^\\c{aQ.), (4.61) 

which in combination with the second equality in fl4.60p then implies ( I4.57P 
by direct verification. 



5. Spherical cloaking device with uniform cloaked contents and 
sharpness of our estimates 

In this section, we consider our near-cloaking scheme within spherical 
geometry and uniform cloaked contents. For this special case, we shall assess 
the cloaking performance, namely Theorem 12. 11 and the result illustrates the 
sharpness of our estimate in Section [31 

In the rest of this section, we choose VL to be Br, R > 0, and a'^ to be a 
scalar constant multiple of the identity matrix, and to be a positive con- 
stant. By a bit abusing of notation, we shall regard cr^ as a scalar constant. 
In the following, we first consider the two-dimensional case. By transfor- 
mation acoustics, it is straightforward to show that 0"^ = cr'^,<la = P'^ 
Dp/2. Let uJa = uj^^Qa/aa and ui = uj^ qijoi = cja/I + ip^^^2 (we choose the 
branch of y/1 + i such that 3\/l + i > 0, that is \/l + i = 2^e*t). Suppose 

oo 
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and according to our earlier discussion in Section 4, 



n=—oo 



(5.1) 



We assume that the solution of f l3.2p is given by 

oo 

{LOa\x\)e'''\ X e Bp/2, 



n=—oo 
oo 



Up{x) = < 



J2 CnJn{ooi\x\)e'^' + J] dnH^'\uJi\x\)e'^' , xeBp\Bp/2, 

.=—00 n=— 00 

00 00 

^ a„J„(a;|x|)e*"V 6„i^(l)(wkl)e'"^ x e BrXB^. 



(5.2) 

We shall denote Ua = Upls^^^^ ui = Up\Bp\Bp^2 = ^\br\Bp- By the 

standard transmission conditions on dBp/2, dBp and the boundary condition 
on dBj^, we have 



diy{x) diy{x) 
ui[x) = ur{x), cri-^-j-^ = ^7773-, X E dBp, 



dunix) 
du{x) 



dv{x) dv{x) 



(5.3) 



TpiyX), X G OBr. 



Plugging the series representations (15. 2p into (15. Sp . we have the following 
linear system of equations for the coefficients, 



'e„J„(wap/2) = c„J„(u;;p/2) + dnH^^\uJip/2), 

(1)' 

^aaqae„J'n{ojap/2) = ^/aiqi[cnJ'n{ojip/2) + d^m {ujip/2)], 
CnJn{(^ip) + dnH^\uJip) = a„J„(u;p) + b„H^\ujp), 
y/aiqi[cnJn{(^ip) + dnH\ {uip)] = anJ'n{up) + {up), 

(1)' 



(5.4) 



Letting A = = ^^ P"^"^' ^^"^ ^^'^ ^^^^ equations of ((53]) we 
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have 



Cn if Jn{uJaP/2) = 0, 



J„{uJip/2) 
H'it\u,ip/2) 

JlMpn)~AMuJlp/2) 2Z''jM 



if J„Kp/2) ^ 0. 



(5.5) 



Denoting the expressions before c„ in (15. 5p by T„, namely (i„ := T„c„, and 
substituting (i„ into the third and fourth equations of (15. 4p . we have by 
straightforward calculations 



/ 4 e 8 p 



1+1 J4(a;;p)+T„Hi^''(a;;p) 



Jn (t^i p) + Tn hJP^UIi p) 



Jniujp) - j;(wp) 



Jn{l^lp)+lnHf^ (l^lp) 



(5.6) 



Let r„ denote the expression before a„ in (15. 6p . namely 6^ := r„a„, and 



^n(p) 



Jn{uJlp) + TnHi^\(^lp) 

Plugging (15.61) into the last equation in (15. 2p . we have 



ijn[Jn{0oR)+TM'\uR)] 



AnO 



whereas the "free space" solution uq{x) E H^{Q) of (I2.15P is 



X e dBji, 



1pnJn{(^\x\) ^ 

ujJ'^{(jjR) 



,in9 



X G Br. 



(5.7) 



(5.8) 



(5.9) 



Hence, 



[up{x) - uo{x) 



oo 



1 + 



T.aH'it^\ujR) 
J^{u,R) 



And 



(5.10) 
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and therefore 



oo ^ 

n — — rv^ ^ 



2\ 1/2 



< 



< 



. n=— oo 

oo 



n=— oo 
1/2 



hr,V2nR 



J v^n=— oo ^ ^ 



uJ'^{uR) 



hr 



Jn{u}R) 



UJ'^UJR) 



H-^/\dBn) 



(5.11) 



where 



1 + 



Since H^_^l{ujr) = {—l)^Hn\ujr), we only need consider n > in estimat- 
ing the series in the last inequality in (15. lip . Using the asymptotic behaviors 

of Jn{z), hI^\z), J'^{z), H!a^\z) as both '^z and ^z tend to +oo (cf. (J, 
(2^). one can show 



(5.12) 



which together with the asymptotic behaviors of Jn{ujp)^ Hn\ujp), J'^{uip), 
H^n^\up) as p ^ +0 (cf. 0,1291), 

one can further show 



p 7TUJI 



[246 8 p2 



f], ^ = 0, 



7ri(a;p)^"[23e*^tjp"+t + r2]/(2"n!)^ n > 1. 



(5.13) 



Then using the estimates in f l5.13p . together with the use of the asymptotic 
developments of the Bessel and Hankel functions for large n (cf. one can 
verify that there exists a sufficiently large integer A^i such that 



p nun 



[2h'fp 



M<-jR) 



Hp''' i^R) -\ ri — f\ 



h ~ irifr :n\'^"^ J^n\^R) _ ''(^fi) l ^ < n < N, 



(5.14) 



K. ~ f (|)^"f2te^f 



e 8 ujp^2 + n\, n > 
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Hence from fl5.14l) . we readily see that there exists a constant Ci independent 
of p for p sufficiently small such that 



\hn\ < Cip^, n < iVi, 



and for n > Ni 



\hJ < 



P_ 
n \R 



2n 



(5.15) 



(5.16) 



Here it is emphasized that due to the asymptotic developments of Hq and 
hi, ( 15.15P is the best estimate one could achieve, namely Cip^ could not be 
improved. Now, using ( ISTTSl) . we see that 



Ni 



n=0 



i?2 ; 



ujJ'„{uR) 



Let A''i be sufficiently large such that 
p < min{i?/4, 1} 



JnjujR) I 



< C2p\ (5.17) 

< 1 for n > Ni, then for 



r 



JnituR) 



< 



/ 

n>Ni ^ ' 

-^y ('i + - 



hr 



n>N-i 



p \ 2(n-l) 

n Kr) 



— 94- — 



(5.18) 



<Csp' 



Combining f l5.1ip . fl5.17p and ( I5.18p . we have 

II^P - ^o||//i/2(aB^) < C'plV^||//-i/2(aB^). 



(5.19) 



Moreover, from the optimality of the estimate (I5.15p . we readily see the 
sharpness of fl5.19p . 

Next, we shall investigate the asymptotic behavior of the boundary con- 



dition on dB'^, namely, -^|a_B . Since 



du 



R 



du 



dBp 



AnO 



where 



In := ttnJ'niuJp) + bnH^^'{up) 
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and from the last equation of fl5.4p 



(5.20) 



uj[J'^{ujR) + VnHl^^ {uR)] 

Hence, we only need study the asymptotic behavior of /„. By direct calcula- 
tions, we have 



2ie*fpi+iH„(p)^„ 



uj[J'^{uR) + T^H^rt^ {uR)] [2ie'fp^+-2nn{p)H^^\ujp) - hI^^ [up]] 

Using the asymptotic behavior of l-Lnip) in (I5.12p . r„ in fl5.13p . and the 
Wronskian J„(t)F^(t) — J'^{t)Yn{t) = ^ (cf. |ll!]), one can show that there 
exists a sufficiently large integer N2 such that 



' 7 i— 
Iq ~ —246 8 p 



1+^ 



^0 



n = 0, 



1 : StT -1 I 5 



0{p 



,71+1+5/2^ 



, l<n<N2, 



(5.21) 



(|)«, n>iV2. 



Using ( ]5.2ip and a similar argument to that for the proof of Lemma 14.21 we 
can show that 



dv 



ip 



(5.22) 



where C is independent of p, 6 and ip. Hence, we readily see that as 5 — +00, 
the lossy layer {Bp\Bp/2;<7i,qi} converges to a sound-hard lajei; that is, the 
normal velocity of the wave filed would vanish on the exterior of the layer, 
namely 



du. 



dv 



ip ■) 



-> as (5 ^ +00. (5.23) 

On the other hand, the sound-hard layer lining is considered in |24], and it 
is shown that one could achieve optimally within the ideal cloaking for 
the regularized cloaking construction. Hence, fl5.22p and f l5.23p also partly 
illustrate the sharpness of our estimates. 

The three-dimensional case could be treated similarly, which we only 
sketch in the following. Let 

00 n 
n=0 m=—n 
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with 

\l'^'H3B.) = E E (1 + + ^)lR'r"'\i^nR? < +00. (5.24) 



00 n 

\ II 9 > 

x) 

n=0 m=—n 



Noting (Ja = cr'a/P: la = Qa/ similar to (15 ■2p for the 2D case, the 
wave fields in the separated domains could be represented as follows 



00 n 

M^)=J2 E On(c^aixi)rr(*), 

71=0 m=—n 

00 n 00 n 



M^)=Y. E cjMx\)y:-{x) + y, E oi'H^d^i)>r(^^), (5.25) 

n=0 m=—n n=0 m=~n 

00 n 00 n 

««(^)=E E cjn(^ixi)rr(*) + E E Khi^\uj\x\)Y:r{s^). 



n=0 m=—n n=Q m=~n 



Similar to (15.31) . using the standard transmission conditions and the boundary 
condition, one could derive the following linear system of equations for the 
coefficients 

f e™j„(a;,p/2) = Cj^Mp^) + d^h'^^\ujip/2) 

CUuJip) + d^h^^\uip) = a^Uoop) + b^hll\iop) ^ (5.26) 
^i[Cf„{uip) + d^h^^^' {uip)] = a^j'ni^p) + C/^i'^Vp), 



Letting A = = V°|^ p~i~^ and solving fl3^ . one has 



where 



j„{uJip/2) 



if jn{0JaP/2) = 0, 



if j„(u;aP/2) ^ 0. 



jn{i^aP/2) 
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Let r„ denote the expression before in f l5.27p . Then 

r 



oo n 



[Up{x)-Ua{x)\\dBR = 5Z 5Z 



Let 



Then 



n=Q m=—n 



h^n\ujR) h'n"{ujR) 



jniulR) jU^R) 



1 + 



Yrix). 



(5.29) 
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1 + 



I / ^ / Ml2 \ - \ - /i , ^(^ + 1) 



.«) = E E + 



?i=0 m=—n 



i?2 



i^l^UujR) 



9nR 



oo n 



V-t I n(n+l) 
H B? — 



. n=0 



n{n + 1) , 

^2 



jn{ujR)gn 



(5.30) 



By similar asymptotic analyses to the 2D case, one can show that there exists 
a sufficiently large integer N^, such that 



go tp [Z-ie s p2UJ —li^^cjBT j,',(<^i?) J' - u, 

^ in+l)(2n)\(2n+l)\\- i„(ujB) j^C^iJ) J' J- ^ /i ^ iV3, 



'(n+l)(2n)!(2n+l)! L j„(w/J) 
p \2n+l 2n+l 



(5.31) 



By a similar argument to the 2D case, applying f l5.3ip to the estimation of 
f l5.30p . one can show that 



\Up{x) -Uoix)\\Hi/2(^QBj,) < Cp^\\^p\\H-l/2(^9BJ,), 



(5.32) 



and moreover, the estimate is optimal. Furthermore, one could also show in 
this 3D case 



du 



R 



du 



//-i/2(aBi) 



< Cp 



1+1 



H-^/^idBR) 



(5.33) 



That is, we also have that the lossy layer would converge to a sound-hard 
layer in the limiting case as 5 — ?■ +00. 
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6. Discussion 



In this work, we consider a novel near-cloaking scheme by employing a 
well-designed lossy layer between the cloaked region and the cloaking reg ion. 



The study follows the spirit of the one developed in [19| . However, in [19| the 



authors rely on a lossy layer with a large lossy parameter for the successful 
near-cloaking construction, whereas we rely on a lossy layer with a large 
density parameter. They are of different physical and mathematical nature. 
As was discussed earlier in Introduction, the lossy layer proposed in 



19 



IS a 

finite realization of a sound-soft layer, whereas the FSH layer in the current 
work is a finite realization of a sound-hard layer. This is confirmed by f l5.22p 
and fl5.33p derived in Section 5, which indicates that as 5 — )■ +oo the FSH 
layer converges to a sound-hard layer. Moreover, we have the following result 
which further supports our such observation. 

Theorem 6.1. Suppose — is not an eigenvalue of the Laplacian on Q\D 
with Neumann boundary condition. Let Ugh ^ H^{VL\D) he the unique solu- 
tion of 



y — 

^ dxj 



N 
N 



du 



X) 



sh 



dx-i 



+ uj'^qp{x)ush 



zn Q\D, 



:du 



sh 



a 



dxj 

dush 
dxi 



^ e H-^/\dn) on on, 



(6.1) 



on dD, 



i)i=l 



and V 



where {fl]a^,q^} is the medium in h2.12^) and, v 
are the exterior unit normals to dD and dVt, respectively. That is, Ush is 
the solution corresponding to a sound-hard obstacle D buried in the medium 
(y^,q^}. Let u G H^{VL) be the solution corresponding to the cloaking 
device, namely. 



N 



N 

E 

k i,j=i 



E ^ ( ^'^(^)'Er- ) + = in 



du 
'dxj 



Via 



du 
dxi 



ip E H-^/^{dn) on dn. 



(6.2) 



36 



Then for sufficiently small p > 0, we have 

\\ush - M||//i/2{9n) < C'P^II^IIh-i/2(9q)> (6-3) 
where C is a constant independent of p and if). 
Proof. Let 

Ush = F*Ush and u = F*u, 

then by transformation acoustics, it is readily seen that u is exactly Up in 
and Ush satisfies 

,2- 



Aush + UJ Ush = in ^\Dp, 
— - — = ip on oil, 

OUsh 



(6.4) 



dv 



on dDr, 



Moreover, we know 



Ush = Ush 



u^h and u = u on dVL. 



Let Uq be the solution of fl2.15p . By straightforward verification, we first 
see that Q = Uq — Ugh £ H^(S^\Dp) satisfies 



Ag + w^g^O in Vl\Dp, 

—— = on oil, 
dQ _ duo 
du du 
By Lemma [4.11 we have 

\\ush - uo\\Hi/2(dn) = II <5 11^1/2(00) < Cp 
On the other hand, by Theorem 13.11 

\\u~uo\\Hi/2(^dn) < Cp^ 
Hence, by (16. 6 p and ( 16. 7p we have 



(6.5) 



N 



H-i/2(an)- 



(6.6) 
(6.7) 



\\Ush — U\\fjl/2(^QQ^ = \\Ush — U\\fJl/2(^QQ^ 

<\\ush — '"o||/fi/2(af7) + \\u — Mo||//i/2(gQ) 
l//-i/2{ar2)- 
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Theorem 16.11 indicates that for a small p, the FSH layer with a large 
density parameter really behaves like a sound-hard layer due to that the 
exterior wave effects are close to each other. Our near-cloaking scheme by 
employing such FSH lining is shown to produce significantly enhanced cloak- 
ing performances compared to the existing ones in literature. The cloaking 
construction is assessed within general geometry and arbitrary cloaked con- 
tents. The assessment is based on controlling the conormal derivative of the 
wave field on the exterior boundary of the FSH layer and estimating the ex- 
terior boundary effects of sound-hard-like small inclusions. Finally, we would 
like to remark that our present study could be extended to the near-cloaking 
of full Maxwell's equations by using the technique developed in this work and 
the estimates due to small electromagnetic inclusions derived in [qJ, which 
will be reported in a future paper. 
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